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NUIJ TYPE PENCILS OF HYPERBOLIC POLYNOMIALS 


KRZYSZTOF KURDYKA AND LAURENTIU PAUNESCU 


Abstract. Nuij’s theorem states that if a polynomial p G M[z] is hyperbolic (i.e. has only 
real roots) then p + sp' is also hyperbolic for any s G K. We study other perturbations 
of hyperbolic polynomials of the form Paiz, s) := p{z) + J2k=i akS^P^^^z). We give a full 
characterization of those a = (oi,..., ad) G for which Pa{z, s) is a pencil of hyperbolic 
polynomials. We give also a full characterization of those a = (oi,..., ad) G R'^ for which 
the associated families pa{z,s) admit universal determinantal representations. In fact we 
show that all these sequences come from special symmetric Toeplitz matrices. 


1. Introduction 


Hyperbolic polynomials and the problem of their determinantal representations is a very 
active area of real algebraic geometry. A nice survey of Vinnikov [9] is a good source on 
recent developments in this subject. The goal of this paper is a study of 1-parameter families 
of hyperbolic polynomials and their universal determinantal representations. Recall that a 
polynomial p G M[z] is called hyperbolic if all its roots are real. Clearly any monic hyperbolic 
polynomial of degree d is a characteristic polynomial of a symmetric dx d matrix. First we 
recall the following theorem proved by W. Nuij [B]. 


Theorem 1.1. Let p G M[ 2 ;] he a hyperbolic polynomial, then 

p -\- sp' 


is hyperbolic for any s G M. 


We give below a proof of this result, based on the existence of determinantal representation 
of the family of the polynomials p + sp' , s G M. In fact we state and prove a generalization 
of Nuij’s result. To this end we propose the following dehnition. 

Definition 1.2. We say that a = (oi,..., a^) G M'’* is a Nuij sequence if for any hyperbolic 
polynomial p of degree d, the polynomial 

d 

Pa{z,s) := p{z) + ^ (1-1) 

k=l 

is hyperbolic for any s G M. We denote by Md the set of all Nuij sequences in M'’*. 
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Note that by Theorem 11.11 a = (1,0,..., 0) is a Nuij sequence for any d & N, d > 1. 
On the other hand, repeated application of Theorem 11.11 also produces Nuij sequences; for 
instance we have 


p + sp' + s{p + sp'y = p + 2sp' + s'^p". 

Hence (2,1, 0,..., 0) is a Nuij sequence for any d G N, d > 2. In Section [3] we shall see 
however that there is an essential difference between those two families, with respect to their 
determinantal representations. 

Surprisingly the set Afd has a nice explicit description. 


Theorem A. A sequence a = {ai,..., ad) E is a Nuij sequence if and only if the polyno¬ 
mial 


qa{z) ■.= z^ + J2 


Ukiz 


)(k) ^ 


k=l 


k=l 


d! 




^d—k 


{d-k)\ 


( 1 . 2 ) 


is hyperbolic. 


In other words, the theorem states that to check that a given a = (oi,..., Od) G is 
a Nuij sequence it is enough to check hyperbolicity of Pa^z., s) only for p{z) = z'^. The 
proof is given in Section [21 it is based on a deep result of Borcea and Branden [T] which 
gives a characterization of linear maps (on the space of polynomials) preserving hyperbolic 
polynomials. 

The second part, developed in Section [31 concerns universal determinantal representation 
of some Nuij sequences. 


Definition 1.3. We say that a = (ai,..., ad) G Afd C admits a universal determinantal 
representation if there exists a symmetric matrix Aa such that for any hyperbolic polynomial 
p of degree d we have 

Pa{z,s) = det{zl D + sAa), (1.3) 

where D is a diagonal matrix whose characteristic polynomial is equal to p = Pa{z, 0). The 
matrix Aa will be referred as a matrix associated to the sequence a = {ai,..., Od). We 
denote by UAfd the set of all Nuij sequences in which admit universal determinantal 
represent at ions. 


Recall that a square matrix is Toeplitz if all parallels to the principal diagonal are constant. 
We say that a symmetric Toeplitz matrix is special if all entries outside the principal diagonal 
are equal to some /5 G M, and of course all entries on the principal diagonal are equal to some 
a G M. In the sequel we will denote such a d x d matrix by Ta^yid), and its determinant by 
ta,/ 3 {d) := detTa^i 3 {d) = {a — (3)'^~^{a + (d — l)/9). We obtain the following characterization 
of all Nuij sequences which admit universal determinantal representations. 

Theorem B. A sequence a = (oi, ..., ad) is a Nuij sequence with a universal determi¬ 
nantal representation if and only if there exit a, /? G M such that 

cij zrta^yii)i 1, . . ., d. 

Z • 

2 


(1.4) 





2. Hyperbolic polynomials and Nuij sequences 

First we recall some facts about the space l-Lf of hyperbolic (monic) polynomials of some 
hxed degree d. For x = (xi,..., Xd) G we have the A:-th elementary symmetric polynomial 

Ck{x)= ^ Xi^---Xi^, ( 2 . 1 ) 

for k = We will identify any b = {bi,... ,bd) G with a monic polynomial 

hb := + Ylt=i Thus we can write Hf = c(M”), where c = (ci ,... ,Cd) : —)• is 

the Viete map, hence by the Tarski-Seidenberg theorem it follows that is semialgebraic. 
Moreover, the Viete map c = (ci,...,Crf) : —)■ is generically a submersion, hence 

Hf = c(M") has nonempty interior. In fact Tif is a basic semialgebraic set which can be 
described using generalized discriminants or Bezoutians (see a nice exposition in [7] or even 
a more detailed in [S]). Recent developments on hyperbolic univariate polynomials are given 
by Rostov in his survey |3]. 

For the proof of Theorem A we need to recall several dehnitions and results from |T]. 

Definition 2.1. [H Dehnition 1] We say that a polynomial f{zi,..., Zn) E C[zi,..., Zn] is 
stable if f{zi ,..., 7 ^ 0 for all n-tuples ( 2 : 1 ,... ,Zn) E C"' with im{zj) > 0, for j = 1,..., n. 

If in addition / has real coefficients, it will be referred to as real stable. The set of stable 
and real stable polynomials in n variables will be denoted by T-LniC) and T-LniM) respectively. 
Note that for n = 1 a polynomial / is real stable, precisely means that / is hyperbolic. 

Let T : Crf[z] —)■ €,d[z\ be a linear map, where Cd[z] stands for the vector space (over C) 
of complex polynomials of degree at most d. We extend it to a linear map T : Cd[z,w] —)• 
Cd[z,w], by setting T{z^w’‘) := T{z^)w^ for all k = 1,... ,d and I E N. We now state the 
result which is crucial for the proof of Theorem [1] 

Theorem 2.2. [H Theorem 4] Let T : Cd[z\ -E Cd[z\ be a linear map. Then T preserves 
stability if an only if either 

(i) T has range of dimension at most one and is of the form 

T{f) = a{f)P, 

where a : CrfRl -E- C is a linear functional and P E 'Hi(C); or 

(ii) T{{z + wy) En2{C). 

Proof of Theorem A. 

Assume a = (ai,..., a^) G is a Nuij sequence. Hence by Dehnition 11.21 applied to 
p{z) = z^ with s = 1 we obtain that the polynomial Pa dehned by fll.ip is hyperbolic. 

To prove the converse let us hx some a = {ai,..., a^) G and assume that the polynomial 
Qa dehned by fll. 2 p is hyperbolic. We associate to the sequence a = {ai,..., Ud) a linear 
operator Ta : Lld[z\ —t Cd[z] dehned by 

d 

Ta{p){z) := p{z) + ^ akP^’'\z) E M[ 2 ;] (2.2) 

k=l 

Lemma 2.3. Ta{{z + w^) = qa{z + w). 
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Proof. We expand first the right-hand side of fl2.2l) 


TMz + wY) ^ T [Y, 


d / ,\ \ ^ / J\ 

j fj—j \ ( ^ \ (j—imr i\ 


.1=0 


i=0 




Note that 


So, 




j=0 


j=0 


= Z^«i 

j=0 


(i-JV- 




i=0 






X */ \' (* — ?)! 

i=o ^ ^ \i=o ^ 


hence 


dl 


%\ 

i=0 ^ ’ \i=o ^ ’ 






On the other hand 


d 


d\ 


\d—i 


+ n;) = —-—ai(z + 

i =0 ^ '' 

The coefficient in (12.31) which comes with a^, j = 0,1,..., d is eqnal to 


(2.3) 


(2.4) 


d\ 


_ — _ 

(d-PU\(i- P\ 


d\ 


^k^d-j-k 


^ {d- i)\i\ {t - j)\ (d-k- j)\k\ 

The coefficient in (12.41) which comes with a^, j = 0,1,..., d is eqnal to 


d! 


{z + wY ^ = 


d! 

i.d-3)\ 


E 

fc =0 


d- j 
k 




E 

i—j=k=0 


d! 


id-k-j)l 




Hence these coefficients are eqnal which proves the lemma. 

□ 


By the assnmption qa has only real roots. Hence qa{z -|- tn) is a stable polynomial in 
variables (z, w). Indeed, if im{z) > 0 and im{w) > 0 then im[z - 1 - tc) > 0, so qa{z -|- tc) 7 ^ 0. 
By Lemma [23] we have Ta{{z -f wY) = daY + w). Applying Theorem 12.21 we conclnde that 
the operator Ta preserves stability, hence Ta restricted to preserves hyperbolicity. Tims 
we have proved that 

d 

Pa{z,l) =p{z) + '^akP^’'\z) 
k=l 
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is hyperbolic whenever p G Mrf[z] is hyperbolic. Let us take s G M* and denote a{s) := 
(sai,..., , s'^Orf). Then the polynomial 

d d I 

qais){z) ■=zUYl + E 

k=l k=l ^ '' 

is again hyperbolic since qa{s){z) = s~‘^qa{sz). Thus, by applying the above argument to 
the sequence a(s), we conclude that Pa{z, s) := p{z) + Yl'k=i <ZkS^P^^\z) is hyperbolic for all 
s G M and any p G hyperbolic. This ends the proof of Theorem A. 

2.1. Iterations of Nuij’s sequences. Let a = (ai,..., a^) G and b = (6i,..., 6^) G 
be two Nuij sequences, we dehne their composition b o a := c = (ci,... ,Cd) in the following 
way. For any polynomial p(z) G M[; 2 ] 


Pc(z, s) = (Pa)b(z, s) = Pa{z, s) + ^ = P + 


Note that with the convention ao = feo = 1 we have 




Let a},... G We dehne by induction the composition of r copies of sequences: 


/i(a^) = a^, Ir{o^ 1 ..., a’’) := ..., a'’ o a*". 

Explicitly, if A(a^,..., a^') = c = (ci,..., q), then 

Ck= (2-6) 

il<-"<irdl-\ - \-ir = k 

Let us consider the original Nuij sequences of the form 

a' = (a:i,0,...,0) gM'", (2.7) 

where Xi G M, i = 1,..., d. Then, /^(a^,..., a*^) = c = (ci,..., q) is the Nuij sequence 
obtained by the iteration of a* and 

Ck= ^ (2.8) 

h<---<ik 

for k = 1,..., d. Thus = Ck{xi,..., Xd) is in fact the k-th. elementary symmetric poly¬ 
nomial of xi,... ,Xd. Denote by c = {ci,... ,Cd) : —)■ the Viete map and recall that 

Pif = c(M”). Thus we obtain that "Hf C J\fd. For d G N let us denote by —)■ the 

following linear map; 


bdi^(^l: • • • 1 ■ ■ ■ ) ®rf) • (d(il, . . . , 


{d-k)\ 


O'k: ■ ■ ■ ) d\cid). 


Theorem A and the above discussion can be summarized as follows. 
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Corollary 2.4. For any d Efi we have 

Example 2.5. For d = 2 we have FLl = {aj — 4 a 2 > 0} C = Wi — 2 a 2 > 0}. 

3. Universal determinantal representations 
We shall consider 1-parameter families of hyperbolic polynomials. A polynomial 

p{z, s) = + ai{s)z'^~^ H-h ad{s) 

will be called a pencil of hyperbolic polynomials if and only if, 

• for each s G M the polynomial z i—)■ p{s, z) is hyperbolic, 

• each coefficient aj(s) G M[s] is of degree at most i. 

For any d > 1, we shall denote by VFLd the space of such pencils of hyperbolic polynomials. 

We say that a polynomial p{z, s) admits a determinantal representation if there are real 
symmetric matrices Ao,Ai such that 

p{z, s) = det(z/ Aq sAi), 

and clearly in this case p{z, s) is a pencil of hyperbolic polynomials. 

As an easy reformulation of a remarkable theorem of Helton and Vinnikov [2] reads 

Theorem 3.1. Any polynomial p{z, s) G VFLd admits a determinantal representation. 

Indeed let us set 2 : = x~^ and s = x~^y and finally 

/(x, y) := x^p{z, s) = x'^p{x~^, x~^y). 

Then / is a real zero polynomial in the sense of Helton-Vinnikov, so it has a determinantal 
representation according to Theorem 2.2 in [2]. In fact, as noticed by Lewis, Parrilo and 
Ramana [5] , Theorem 13.11 is a positive answer to the nonhomogeneous version of the Lax 
conjecture [ 1 ]. 

We want to characterize all Nuij sequences a = (oi,..., a^) G such that for any p G M[;2], 

hyperbolic polynomial of degree d, the associated pencil of hyperbolic polynomials 

d 

Pa{z,s) :=p + y^afcgV^^ G R[z], 
k=l 

admits a universal determinantal representation] by this we mean that there exists a sym¬ 
metric matrix Aa such that 

Pa{z, s) = det{zl + D + sAa), 

where D is a diagonal matrix. In other words —D has on the diagonal all the roots of p 
written in an arbitrary order. The matrix Aa will be referred as a matrix associated to the 
sequence a = (oi,... , 0 ^). We denote by UAfd the set of all Nuij sequences in which 
admit universal determinantal representations. 
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3.1. Special Toeplitz matrices. Recall that a square matrix is called Toeplitz matrix if 
all parallels to the principal diagonal are constant. We say that a symmetric Toeplitz matrix 
is special if all entries outside the principal diagonal are equal to some G M, and of course 
all entries on the principal diagonal are equal to some a G M. We will denote such a matrix 
by Ta,/3. 

In the next proposition we will show that special Toeplitz matrices give all Nuij sequences 
which admit universal determinantal representations. 

Proposition 3.2. Let a = (oi,..., a^) G UMd- Then, there exists a special Toeplitz matrix 
Ta,i 3 which is associated to the sequence a. The constant a is unique. For d = 2 we have two 
choices [3 or —/3. If d> 3, then /? is uniquely determined. 

Proof. Let us £x a sequence a = (ai, ... ,ad) G UMd, and let Aa be a symmetric matrix 
associated to a. It means that for any hyperbolic polynomial p G M[z] we have 

Pa{z,s) = d(it{zl + D F sAa), (3.1) 

where D is a diagonal matrix with characteristic polynomial equal to p. We will hnd a 
special Toeplitz matrix T^^p such that 

Pa{z, s) = det( 2 ;J + D + sTa^y). 

As a piece of convention, we recall that a, j x j minor of Aa is principal if it is the 
determinant of a matrix obtained from Aa by deleting rows and columns containing d — j 
elements from the principal diagonal. With the assumption of Proposition 13.21 we have. 


Lemma 3.3. For any j = 1,... ,d, all j x j principal minors of Aa are equal. 


Let — Ai,..., —Xd be the roots of p. Since p can be chosen arbitrarily, we may consider 
both sides of the identity l3.ll as polynomials with real coefficients in variables Wi \= z + Aj, 
i = l,...,d. Since R is a held of characteristic 0, the coefficients corresponding to the 
monomials in Wi^ ■ ■ ■ Wi., where R < • • • < fj, on right-hand and left-hand side are equal. It 
is enough to expand both sides to check the statement of the lemma. In particular the 1x1 
minors, which are actually the entries on the principal diagonal, are all equal to some a G R. 


Lemma 3.4. Let Aa = {aij). Then there exists /3 G R such that for any distinct i,j we have 
a% = fiT 


Indeed to each entry a^, i ^ j we can associate the 2x2 principal minor 


det 


a 

Qij 



Hence by Lemma 1X51 all oA are equal for i ^ j. We put (3"^ = al. Clearly the statement 
of Proposition 13.21 is trivial for /3 = 0, so in the sequel we assume that (3^0. 

Before analyzing the case of j x j principal minors, where j > 3, we need an explicit 
formula for the determinant of a special Toeplitz matrix Ta^y. 


Lemma 3.5. IfTa^y is a special Toeplitz matrix of size d x d, then 

taAd) ■= det Ta^p = {a- I3f~^{a + {d - l)/3). 
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(3.2) 






Next we consider the 3x3 principal minors of the matrix Aa- We know by Lemma f3.4l 
that for any i ^ j we have aij = eij\l3\, where eij G {—1,1}. We will show that the sign of 
eij can be nniformly chosen, which means that either eij = 1 for all i ^ j, or etj = — 1 for all 
i ^ j- Le ns write this minor in the form 


/ a epl/5| eik\l3\ \ 

det I eij\^\ a ejk\/3\ \ = + 2eijeikejkP^\P\ - 3al3^. (3.3) 

\ (^ik\l3\ ejk\l3\ a j 

By Lemma 13.31 all these minors are eqnal, so there exists ^ G { — 1,1} snch that for all 
choices \ <i < j <k < dwe have 

( 2 '^) 

This shows that we may chose e^j = ^ for all i ^ j- 

Assnme now that d > 4. We have to show that if we pnt e^- = ^ for any i ^ j, then 
actnally all principal minors j x j, j > 4: are eqnal to the valne of a principal minor j x j, 
j > 4 for the original matrix Aa, so in fact they are determined jnst by Note that it is 
enongh to consider the case a = 0 and /3 = 1. First we consider the case d = 4, so 

ei3 ^14 ^ 

^23 ^24 

0 634 

£34 0 J 

by eii and nse relation 13.41 Thns we obtain 

^13 ^14 ^ 

■^^13 ^^14 

0 ^ei4 

^ei3 0 / 

)f Ba by eij and nse the fact that = 1. So 

1 1 \ 

0 e ■ 

^ 0 

Multiplying the hrst row and the hrst colnmn of Ca by we can see that 

det Ca = det To,i = to,i(4) = -3. 

Bnt on the other hand det Ca = ( 612613614 )^ det Aa = det Aq. Accordingly we may assnme 
that Aa = The same argnment applies for any d > 4. Hence the existence in Proposition 
13.21 follows. To proof the nniqneness, note that a and are nniqnely determined. Clearly 
the eqnation 03 = 4 ^( 0 ; — /5)^(a + 2/5) nniqnely determines /5. □ 

As a conseqnence we obtain the following characterization of Nnij seqnences which admit 
nniversal determinantal representations. 


Aa = 


/ 


V 


0 

ei2 

ei3 

ei4 


ei2 

0 

^23 

^24 


For each i > 2 we mnltiply the Ah row of Aq 
the matrix 


Ba-.= 


/ 


V 


^12 

0 

^ei2 

^ei2 


For each j > 2 we mnltiply the jth colnmn c 
we obtain the matrix 

/ 0 1 

Ca := ^ ° 


V 














Theorem B. A sequence a = {ai,, a^) is a Nuij sequence with a universal determi- 
nantal representation if and only if there exits a, /? G M such that 

tti = \t^^y(i), i = 1,... ,d. (3.5) 

i\ 

Proof. If is a special Toeplitz matrix, then for any hyperbolic polynomial p{z) = [z + 
Ai)... (z + Xd) we have a pencil of polynomials 

d 

Pa{z,s) := p + = det{zl + D + sT^^y) 

k=l 

where a* = ^ta,i 3 {i), and D is a. diagonal matrix with entries Ai,..., A^. So the sequence 
a = (oi,..., Od) is a Nuij sequence with a universal determinantal representation. Conversely, 
if a = (tti,... ,ad) G is a Nuij sequence with a universal determinantal representation, 
then by Proposition 13.21 the associated matrix can be chosen as a special Toeplitz matrix 
Hence ai = \ta,f}{i). □ 

Example 3.6. Note that the original Nuij sequence a = (1, 0,..., 0) has a universal determi¬ 
nantal representation. Indeed, Ti^i, which has all entries equal to 1, is the matrix associated 
to this sequence. Note, that this also proves Nuij’s Theorem ll.il 

Remark 3.7. A composition of the original Nuij sequence a = (1,0,.. .,0) with itself gives 
a Nuij sequence b = (2,1, 0,..., 0) which has no universal determinantal representation for 
d > 3. Indeed, if there exist a,/? G M such that R = j]ta,i 3 {i), i = 1,2,3, then a = 2 and 
— 13“^ = 2. Hence /3 = ±\/2. But, then 663 = + 2(3^ — 3al3‘^ 7 ^ 0, so 63 7 ^ 0, which is a 

contradiction. 
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